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Magneto-transport of interacting electrons in a Sierpinski gasket fractal is studied within a mean 
field approach. We work out the three-terminal transport and study the interplay of the mag- 
netic flux threading the planar gasket and the dephasing effect introduced by the third lead. For 
completeness we also provide results of the two-terminal transport in presence of electron-electron 
interaction. It is observed that dephasing definitely reduces the transport, while the magnetic field 
generates a continuum in the transmission spectrum signaling a band of extended eigenstates in 
this non-translationally invariant fractal structure. The Hubbard interaction and the dephasing 
introduced by the third lead play their parts in reducing the average transmission, and opens up 
gaps in the spectrum, but can not destroy the continuum in the spectrum. 
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I. INTRODUCTION 

Electronic transport in low dimensional systems is a 
key to explore important interference effects, local cur- 
rents, switching mechanisms and other unique properties 
that are pre-requisites for estimating the potentiality of 
these systems as quantum interference devices [ll-Q- ^ 
considerable amount of work in this direction has already 
revealed the unique features of phase coherent electron 
transport through quantum dots (QD), Aharonov-Bohm 
(AB) rings, and model molecular systems [3-[l9j. 

A major part of the studies in low-dimensions so far 
describes electron transport in the so called two-terminal 
nano- or mesoscopic devices [l^l- Indeed, this is a fast de- 
veloping field, and has stimulated lot of theoretical work 
based on the non-equilibrium Green's function approach 
within the density functional theory [2l| - |2^ . 

Comparatively speaking, much smaller volume of lit- 
erature on the three or four terminal electronic transport 
have come up in recent times [lll - fisj . A two-terminal 
device is essentially a single path device. The transport 
here is marked by a sharp jump in the transmission phase, 
and is constrained by the Onsager relations of time rever- 
sal symmetry [25j and the current conservation [lo| . So, 
the inclusion of a third terminal that allows the current 
to flow out of the system and breaks the unitarity condi- 
tion, is likely to be useful in extracting useful information 
related to the quantum coherence is low-dimensional sys- 
tems [li,!!^. 

In this communication we undertake an in-depth study 
of the three-terminal magneto-transport of interacting 
electrons in a fractal network. Specifically, we choose a 
fractal geometry following a Sierpinski gasket (SPG) [28l - 
[36j . Such a planar gasket is shown in Figs. [T] and [U 
and can be thought to be equivalent to a self-similar ar- 



rangement of single level QD's [37|,|38[ sitting at the ver- 
tices of each elementary triangle. With the present day 
advancement in lithography practically any design can 
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FIG. f : (Golor online). A 3-rd generation Sierpinski gasket is 
attached to two semi-infinite ID metallic leads, viz, lead-f and 
lead-2. Each elementary plaquette of the gasket is penetrated 
by a magnetic flux 4>- The filled red circles correspond to the 
positions of the atomic sites. 

be tailor made, and the present work thus offers an ex- 
cellent opportunity to study the simultaneous effects of 
magnetic field, electron-electron interaction and the de- 
phasing caused by the introduction of a third lead in the 
system. 

Our motivation behind the present work is two-fold. 
First, we observe that the multi-terminal transport in 
systems with multiply connected geometry is a very little 
addressed (or, unaddressed) problem. In particular, the 
in-built self similarity of systems such as the SPG opens 
up the possibility of investigating the tunneling or switch- 
ing aspect of these systems at arbitrarily small scales of 
energy. Secondly, it is essential, for a completeness in the 
understanding of the spectral properties of fractal net- 
works to know, if the well known multi-fractal. Cantor 
set energy spectrum of non-interacting electrons p8l - [36| 
still retains its character even in the presence of electron- 
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electron interaction or dephasing caused by a third elec- 
trode. In a recent work, the effect of electron-electron 
interactions on the persistent current in a closed loop 
SPG has been reported [s^, but no result exists for open 
self-similar systems connected to electron reservoirs by 
multiple leads. 

On its own merit, the effect of electron-electron inter- 
action on the spectral properties is of great importance. 
Several experiments done on fractal networks have stud- 
ied the magnetoresistance, the superconductor-normal 
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FIG. 2: (Color online). A 3-rd generation Sierpinski gasket is 
attached to three semi-infinite ID metallic leads, viz, lead-1, 
lead-2 and lead-3. Each elementary plaquette of the gasket 
is penetrated by a magnetic fiux 0. The filled red circles 
correspond to the positions of the atomic sites. 

phase boundaries on Sierpinski gasket wire networks (40l - 
l44l |. These experiments pioneered the actual observa- 
tional studies of spectral properties and flux quantization 
effects on planar networks and the Aharonov-Bohm effect 
in systems without translational invariance. Although in 
an early paper the problem of interacting electrons on a 
percolating cluster that displays a fractal geometry , 
has been addressed, to the best of our knowledge, no rig- 
orous effort has been made so far to unravel the effect of 
an interplay of electron-electron interaction and an ex- 
ternal magnetic field on deterministic networks such as a 
Sierpinski gasket (SPG), even at a mean field level. 

Thus, apart from a critical investigation of the possi- 
bility of a fractal device, the present work is also likely 
to throw light on the fundamental spectral properties of 
the deterministically disordered systems. 

We find quite interesting results. Working within a 
tight-binding framework we develop a mean field method 
of studying the three-terminal transport in interact- 
ing systems. The method is then applied to a planar 
SPG and we show that, dephasing definitely reduces the 
corner-to-corner propagation of electrons, with or with- 
out the electron-electron interaction. The magnetic field, 
in the absence of the electron-electron interaction gener- 
ates an apparent continuum in the transmission spec- 
trum. The continuum has already been observed in sev- 



eral fractal lattices [iel - lisj including an SPG [ill, and 
is practically unaffected even when the electron-electron 
interaction is switched 'on'. A detailed study on the ef- 
fect of positioning of the third electrode is also made, 
and the transport properties of an anisotropic gasket are 
compared with its isotropic counterpart. For a compar- 
ative study, the results of the two-terminal transport are 
also presented along with the three-terminal cases. 

In what follows, we present our model quantum system 
in Section II. The essentials of the mean field calculation 
are discussed in Section III. Section IV contains the nu- 
merical results and related discussions, and we draw our 
conclusions in Section V. 



II. THE MODEL QUANTUM SYSTEMS 

Let us begin by referring to Fig. [1] where a 3-rd genera- 
tion SPG is attached to two semi-infinite one-dimensional 
(ID) metallic leads, namely, lead-1 and lead-2 via the 
atomic sites A and B. Each elementary plaquette of the 
gasket is threaded by a magnetic flux (j) (measured in unit 
of the elementary flux quantum = ch/e). The filled 
red circles correspond to the positions of atomic sites in 
the SPG. In a Wannier basis the tight-binding Hamilto- 
nian of the interacting gasket reads, 

^SPG = J2 ««'^cLc»^ + ^ V \f^Aa^3'y + 

+ 1] ^"=lt°'tC4C*4 (1) 

i 

where, is the on-site energy of an electron at the 
site i of spin a (f, \) and v is the nearest-neighbor hop- 
ping strength. In the case of an anisotropic SPG, the 
anisotropy is introduced only in the nearest-neighbor 
hopping integral v which takes on values and Vy for 
hopping along the horizontal and the angular bonds, re- 
spectively. Due to the presence of magnetic fiux (j), a 
phase factor 9 = 27r(/)/3 appears in the Hamiltonian when 
an electron hops from one site to another site. A negative 
sign comes in when the electron hops in the reverse di- 
rection. As the magnetic filed associated with does not 
penetrate any part of the circumference of the elemen- 
tary triangle, we ignore the Zeeman term in the above 
tight-binding Hamiltonian (Eq. [1} . and Cia are the 
creation and annihilation operators, respectively, of an 
electron at the site i with spin a. U \b the strength of 
on-site Coulomb interaction. The Hamiltonian for the 
non-interacting leads can be expressed as, 

^lead = focfc, + (c\cj + c]c,) (2) 

i <ij> 

where different parameters correspond to their usual 
meaning. These leads are directly coupled to the gas- 
ket where the hopping integral between the lead-1 and 
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gasket is ti, and, it is T2 between the gasket and lead- 
2. With this setup we investigate two-terminal electron 
transport through an SPG. 

For three-terminal quantum transport we connect an 
additional lead with the gasket. A schematic view of a 
3-rd generation SPG attached to three semi-infinite fD 
metallic leads, viz, lcad-1, lcad-2 and lead-3 is shown 
in Fig. [21 The gasket and the side-attached leads arc 
described by the same prescriptions as described above. 



III. THE MEAN FIELD APPROACH 



spin Hamiltonians. Then we calculate a new set of values 
of (nil) and {ui^). These steps are repeated until a self 
consistent solution is achieved. 



C. Two-terminal quantum system 

Now we are at the stage of calculating electron con- 
duction across an SPG. 

To determine two-terminal conductance (q) of the gas- 
ket, we use Landauer conductance formula j20l |. At much 
low temperature and bias voltage it can be expressed as. 



A. Decoupling of the interacting Hamiltonian 

Before going to the calculation of electronic trans- 
mission probability through the interacting model of an 
SPG described by the tight-binding Hamiltonian given in 
Eq.[Tl first we decouple the interacting Hamiltonian using 
the generalized Hartree-Fock approach [sol - fs^ . The full 
Hamiltonian is completely decoupled into two parts. One 
is associated with the up-spin electrons, while the other 
is with the down-spin electrons. The on-site potentials 
get modified appropriately, and arc given by. 



(3) 



(4) 



where, n^o. = cl^Cia is the number operator. With these 
site energies, the full Hamiltonian (Eq. [ij can be written 
in the decoupled form (in the mean field (MF) approxi- 
mation) as. 



MF = H '^^t'^'t + " [' 

i 

= ^SPG.t + -^SPG.4 - H U{n,^){n^i 



(5) 



where, -ffgpQ and -ffgPG i correspond to the effective 
tight-binding Hamiltonians for the up and down spin 
electrons, respectively. The last term is a constant term 
which provides a shift in the total energy. 



B. Self consistent procedure 

With these decoupled Hamiltonians (ffgpQ ^ and 
iJgpQ ^ of up and down spin electrons, now we start our 
self consistent procedure considering initial guess values 
of (ni-t-) and {riii). For these initial set of values of {rii-^) 
and {nil), numerically diagonalizc the up and down 
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where, and Tj^ correspond to the transmission proba- 
bilities of up and down spin electrons, respectively, across 
the SPG. Since no spin-flip scattering term exists in the 
Hamiltonian (Eq. [1} , spin-flip transmission probabilities 
will not appear in Eq. [6l In terms of the Green's func- 
tion of the gasket and its coupling to side-attached leads, 
transmission probability can be written in the form [20j , 



T„ - Tr Fi GgpG^^ T2 GgpQ_^ 



(7) 



where, Fi and F2 describe the coupling of the SPG to 
the lead-1 and lead-2, respectively. Here, GgpQ and 

GgpQ ^ are the retarded and advanced Green's func- 
tions, respectively, of the SPG including the effects of 
the leads. Now, for the complete system i.e., the SPG 
and two leads, the Green's function is expressed as. 



Go- — (E — Hrj 



(8) 



where, E is the energy of the injecting electron. Evalu- 
ation of this Green's function needs the inversion of an 
infinite matrix, which is really a difficult task, since the 
full system consists of the finite size gasket and two semi- 
infinite ID leads. However, the full system can be parti- 
tioned into sub-matrices corresponding to the individual 
sub-systems and the Green's function for the gasket can 
be effectively written as. 



G 



SPG. 



^ - ^SPG.a - Si - Si 



(9) 



where, Ei and E2 are the self-energies due to coupling 
of the gasket to the lead-1 and lcad-2, respectively. All 
information of the coupling are included into these self- 
energies. 

D. Three-terminal quantum system 

In order to calculate the conductance in three-terminal 
SPG, we use Biittikcr formalism, an elegant and simple 
way to study electron transport through multi-terminal 
mesoscopic systems. In this formalism we treat all the 
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leads (current and voltage leads) on an equal footing and 
extend the two-terminal linear response formula to get 
the conductance between the terminals, indexed by p and 
q, in the form [13, S^, 



9pq 



(10) 



where, Tpq^^ gives the transmission probability of an elec- 
tron with spin a (f, i) from the lead-p to lead-g. Now, 



> 






-3 3 
E 


- 
- 

, 1 






-3 3 

E 


- 
- 

" ,1 jJ 






FIG. 3: (Color online). Two-terminal conductance g as a, 
function of energy E for a 4-th generation isotropic {vx ~ 
Vy = 1) Sierpinski gasket (A'^ = 42), where lead-1 and lead-2 
are connected at the positions A and B, respectively. The first 
and second columns correspond to <j) = and 0.25, respec- 
tively, where the upper, middle and lower panels represent 
U — 0, 1 and 2, respectively. 

similar to Eq. [7] the transmission probability Tpq^^ can 
be expressed in terms of the SPG-lead coupling matrices 
and the effective Green's function of the SPG as [2(1, 



temperature and use the units where c = h = e = 1. 
Throughout our numerical work we set e^-j- = = for 
all i and choose the nearest-neighbor hopping strength 
V = 1. In the anisotropic case we select Vx = 1 and Vy = 2 
throughout. For side-attached leads, the on-site energy 
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FIG. 4: (Color online). Two-terminal conductance g as a 
function of energy E for a 4-th generation anisotropic {vx ~ 1 
and Vy = 2) Sierpinski gasket (A^ = 42), where lead-1 and 
lead-2 are connected at the positions A and B, respectively. 
The first and second columns correspond to = and 0.25, 
respectively, where the upper, middle and lower panels repre- 
sent U — 0, 1 and 2, respectively. 

(eo) and nearest-neighbor hopping strength (to) are fixed 
at and 1.5, respectively. The hopping integrals (ri, T2 
and Ts) between the leads and SPG arc set at 0.5. Energy 
scale is measured in unit of v. All the essential features 
of electron transport are obtained both for an isotropic 
gasket and its anisotropic counterpart. 



T — Tr 



^ P *^SPG,a ^ 9 "^SPG.a 



(11) 



In the presence of multi-leads, the effective Green's func- 
tion of the SPG becomes (extension of Eq. [9]) [20| . 



'^SPG,^ = \ E-H< 



SPG„ 



(12) 



where, Sp is the self-energy due to coupling of the SPG 
to the lead-p and the sum over p runs from 1 to 3. 

In the present work we inspect all the essential features 
of magnetic response of an SPG network at absolute zero 



IV. NUMERICAL RESULTS AND DISCUSSION 
A. Two-terminal quantum transport 

The two-terminal conductance of a 4-th generation 
isotropic gasket is shown in Fig. |3] with the leads con- 
nected in the positions A and B only. The left panel 
shows the conductance in zero magnetic field, and rep- 
resents the familiar fragmented, scanty distribution of 
the transmission resonances that mark such a lattice. 
The on-site Hubbard interaction displaces the resonance 
peaks, but no marked changes in the spectrum is ob- 
served. 
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With the magnetic field turned 'on' there is how- 
ever, a remarkable change. The spectrum apparently 
exhibits continuous parts, which have previously been 
reported to support extended single particle states 
for spinless, non-interacting electrons. We observe here 
that, the electron-electron interaction, though reduces 
the overall conductance of the system, preserves the con- 
tinuum at the central part of the spectrum. As the in- 
teraction U is increased, the central continuum is more 
or less undisturbed, though the average transmission am- 
plitude still exhibits lower values compared to its J7 = 
counterpart. In addition, there is a signature of opening 
up of new gaps in the spectrum as it appears in Fig.[3jf). 

In Fig. m we show the two-terminal conductance across 
an SPG of the same size as before, but now with a hop- 
ping anisotropy such that, = 1 and Vy = 2. The 
effect of hopping anisotropy is, in general, to increase 
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FIG. 5: (Color online). Three-terminal conductance gi2 
as a function of energy E for a 4-th generation isotropic 
{vx = Vy = 1) Sierpinski gasket (A*' = 42), where lead-1, 
lead-2 and lead-3 are connected at the positions A, B and 
C, respectively. For this arrangement of side attached leads 
(symmetric configuration), all gij's {i ^ j) are identical to 
each other, (a) U ^ 0, (/> = 0; (h) U ^ 0, (f) ^ 0.25; (c) U = 2, 
= and (d) [/ = 2, <;/> = 0.25. 

the overall conductance. The sparse spectrum in the left 
panel of Fig. [3] is now replaced by a substantial density 
of resonant transmission peaks. The introduction of U 
still reduces the average conductance, but the action is 
now somewhat 'delayed'. We need a larger value of U 
to generate a noticeable change in the conductance spec- 
trum. This is evident from the left panel of Fig. |4l On 
the right panel, the interplay of the magnetic field and 
the electron-electron interaction is displayed. The central 
continuum of the isotropic gasket case is still there, but 
with numerous conductance minima somewhat obscuring 
the continuity of this part. As U increases, one observes 



genuine anti-resonances in the previously sustained con- 
tinuum in the spectrum. 



B. Three-terminal quantum transport 

We now discuss the three-terminal transport in a SPG 
network. On the left panel of Fig. [S] the results of the 
conductance between the leads 1 and 2 are shown. The 
third terminal 3 connects to an electron reservoir. The 
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FIG. 6: (Color online). Three-terminal conductance gij as 
a function of energy E for a 4-th generation isotropic {v^ = 
Vy = 1) Sierpinski gasket {N = 42), where lead-1 and lead-2 
are connected at the positions A and B, respectively. Lead- 
3 is attached at the centre of the line BC. This is the so- 
called asymmetric configuration of the leads. First column 
corresponds to (f) — 0, while the second column represents 
(j) = 0.25. For all these spectra we set (7 = 2. 



effect of dephasing is obvious. The spectrum in zero mag- 
netic field is even thinner now, and the interaction U 
plays its part. The magnitude of the resonance peaks 
get suppressed compared to the coherent transport in 
the two-terminal case. A fine scrutiny (not shown here) 
will reveal a level broadening due a dominant phase ran- 
domizing effect [53| leading to a loss of phase coherence. 
Strangely enough, in presence of the magnetic field, the 
dephasing effect is just not enough to destroy or alter 
appreciably, the central continuous part of the conduc- 
tance spectrum. We have carried out extensive numerical 
analysis of the continuum in the conductance spectrum, 
and it appears that the band remains intact over finer 
scales of energy interval and even for a larger sized SPG. 
Therefore, we arc tempted to conclude that there exists 
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a band of extended eigenstates (characterized by finite 
transmittivity) in an SPG fractal threaded by a mag- 
netic field perpendicular to the plane of the gasket and 
that, its a quite robust band, unperturbed by either the 
electron-electron interaction, or the presence of any de- 
phasing effect. 

In the asymmetric case however, the conductance is 
strongly sensitive on the lead positions. For example, 
with Vx = 1 and Vy = 2, the conductance gi2 across the 
1-2 terminals is much larger compared to 323 or gi^ , par- 
ticularly in presence of the magnetic field. The apparent 
continuum in the spectrum persists, as before, even with 
an increasing value of U. The results are displayed in 
Fig. El 

V. CONCLUSION 

In conclusion, we have examined in details, the three- 
terminal transport in a planar Sierpinski gasket fractal 



with a magnetic field threading each elementary trian- 
gle. The electron-electron interaction is considered in 
the Hubbard form, and the Hamiltonian is solved within 
a generalized Hartree-Fock scheme. The sensitivity of 
the conductance spectrum on the positioning of the leads 
is studied extensively, and some prototype results have 
been displayed and discussed. It is seen that the major 
role of the third lead is to bring down the average con- 
ductance of the system. The same role is played by the 
Hubbard interaction. The anisotropic gasket has been 
found to be more conducting than its isotropic counter- 
part with or without the electron-electron interaction. 
However, the curious result is that, the continuum of 
states generated in the spectrum of an SPG as the mag- 
netic field is turned on, is practically unaffected by the 
dephasing effect caused by the additional electrode or by 
the on-site Hubbard interaction. 
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